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Abstract 

We calculate the correlation functions of uniform magnetization in thin ferromag- 
netic films for small deviations from equilibrium, by using a functional formalism. To 
take account of dissipation and fluctuations consistently, the magnetization is coupled 
to a bosonic heat bath. The correlation functions show strong dependence on the na- 
ture of the coupling between the bath and the system. Depending on what coupling we 
choose, we show how the recent results (J. Appl. Phys. 90, 5768(2001); Phys. Rev. B 65, 
172417(2002)) obtained by macroscopic methods can be related to the microscopic 
treatment adopted here. 



1 Introduction 



The problem of magnetic noise in nano-systems and in particular in giant-magnetoresistive 
(GMR) heads is of considerable importance to the physics of magnetic recording. For macro- 
scopic systems, the Landau-Lifshitz-Gilbert equation (LLG) has been used frequently to 
study thermal fluctuations in magnetic systems mainly through a stochastic approach. (Jl|) 
Magnetic noise in thin anisotropic films has been recently treated by at least a couple of 
different methods. (@; H 0) The first work (0) is based on a linearized LLG and application 
of the fluctuation dissipation theorem (FDT).(0) The second method (0) is also approxi- 
mate and based on analogies with the harmonic oscillator. (0) Both calculations are classical, 
however they give different answers. Further, in both works it is assumed that the system is 
close to equilibrium. To better understand this discrepancy, we study this particular prob- 
lem of magnetic noise in thin films from a microscopic point of view. Since we are dealing 
with a magnetization slightly disturbed from equilibrium, we use bosonic degrees of freedom 
to describe the magnetization. (0) Moreover, we use the language of coherent states (CS) to 
describe the states of the magnetization. (|9|) The dissipation is simulated by coupling the 
magnetization to a bosonic environment. (llOl ) The noise spectrum is found by calculating 
the correlation functions of the magnetization. We use a functional method borrowed from 
Field Theory to carry out the general calculation. Eil ) These methods are attractive 
because they are equally applicable to highly non-equilibrium situations and very suitable 
for the many-body problem. Two of us have already used these methods in a recent paper 
that addressed the conditions under which a LLG equation can be recovered from a simple 
quantum model, (flih The results presented here complement those presented in Ref. fjl3h . 
However this paper can be read independently of our previous work. The major result of 
this work is a general expression for the correlation functions from which we can recover 
both LLG-type correlation functions and oscillator-like correlation functions. Currently ex- 
periments do seem to favor the LLG result however we will not address these questions here. 



(JJ; Il4r ) We simply show that for systems close to equilibrium, we can have different theories 
for a large spin. However, it must be kept in mind that the results obtained are based on 
a very simple coupling between the magnetization and the bath. A more realistic coupling 
such as that of conduction electrons interacting with localized magnetic moments is treated 
elsewhere, fflil ) In this latter case the use of functional methods is well justified. 

The paper is organized as follows. In Sec. II, we first introduce the Hamiltonian for the 
system considered here. Then we write this Hamiltonian in an approximate form valid for 
near-equilibrium cases. It is shown that in this case the spin Hamiltonian is exactly that 
of an oscillator. Since correlation functions of the damped harmonic oscillator in a thermal 
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bath are available in the literature, we deduce from them the corresponding correlation func- 
tions for the magnetization. This result does not apply to LLG-type correlation functions, 
however, which are deduced from a different coupling between the spin and the bath. In Sec. 
Ill, we review the harmonic oscillator CS and their relation to those of spin CS. We also 
write the CS bosonic propagator at finite temperature to be used in subsequent sections and 
introduce a generating functional for these propagators in real-time. In Sec. IV, we derive 
the CS generating functional for the system considered here. In Sec. V , a normal mode 
analysis of the noise is carried out. Here we recover the Safonov-Bertam 8 result in the 
limit when the range of frequencies are around the resonance (FMR) frequency. In Sec. VI, 
we derive the general correlation functions for the magnetization. The LLG result is shown 
to follow from the general result by assuming a bath for which the product of the density 
of states and the coup ling constants is linear with frequency. This is the same condition 
recovered in Ref. (|13l ) and is independent of the Hamiltonian of the spin system alone. In 
Sec. VII, we summarize our results. Finally in an appendix we show how these results can 
also be obtained form methods of equilibrium statistical mechanics. 

2 Quantum Model 

In this section, we introduce the quantum model and approximate the magnetization 
operator by a Bosonic operator stressing the analogies with the harmonic oscillator. This 
picture is in general true for close to equilibrium states, irrespective of the equation of motion 
we are using, i.e., LLG or others. 

We consider the following system: A thin magnetic slab with easy axis along the z-axis 
and a hard axis along the x-axis which are in-plane. We assume that there is a large external 
magnetic field H along the easy axis that keeps the average magnetization in-plane. We will 
be interested only in fluctuations around the equilibrium position of the magnetization, i.e., 
fluctuations in the x and y components of the magnetization. The Hamiltonian for such a 
system has the general form (h = 1) 

H = -HS z -K^ + K 2 Sl + J2"kbtb k + v(b k ,blsy (1) 

k 

where K\ and Ki are the anisotropy constants. uj k is the energy of the k-th bath's oscillator 
and 7fc's are the coupling constants, which can be time-dependent. The spin-bath interaction 
V will be taken linear in S and the bath variables b k . A possible form for V is 

V (b k , b[, S)=J2 (ilS+bh + IkbtS-) . (2) 
k 
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The Bosonic-type operators for the bath oscillators, bk satisfy the usual commutation rela- 
tions, 



For the magnetization, S, we have the following commutation relation 



2S Z , 



where 



S + 



S x + iSy, 
: S x iSy. 



(3) 

(4) 

(5) 
(6) 



K\ and K 2 are the anisotropy constants. Uk is the energy of the fc-th oscillator and 7fc's 
are the coupling constants, which can be time-dependent. The Hamiltonian is therefore 
comprised of three terms: the magnetization (spin term), the spin- reservoir interaction and 
the reservoir, 

H = Hs + H-sr + Hr- (7) 
If S is the magnitude of the magnetization, we are interested in the case in which 



S-S 2 



2S 



«1, 



(8) 



that is deviations from the z-axis are small. For precession around the equilibrium position, 
i.e., the z-axis, S z can be assumed constant and the Hamiltonian expression for the spin part 
can be simplified to be of the general form 



where 



2 v ^ ^ + B S y 



(9) 



A = K 1 + K 2 , B = Kl (10) 

On average, we have (S x (t)) = (S y (t)) = 0. To account for thermal fluctuations, we need to 
calculate the two-point correlation functions of the components of the magnetization, i.e., 

(11) 



S x (t)S x (0)),{S y (t)S x (0) 

These are the principal quantities that are needed for a full account of the noise or dissipation 
of the magnetization problem for the model considered here. Since the Hamiltonian is 
quadratic, all higher correlation functions are zero. 



3 



It is now more appropriate to define the following operators, 

1 



and 

then, we have 

If we normalize by 2S, we simply have 



(25) 1 
1 



-S-, 



and 



(25) ! 
S z ~ S — a + a 

a+ = S_ 

5c = - (a + a + ) , 
= — (a - a+) . 



Therefore 



2 



<< 1 . 



(12) 

(13) 
(14) 

(15) 
(16) 

(17) 
(18) 

(19) 



The operators a and a + then behave as bosonic degrees of freedom, i.e., the magnetization 
behaves, in this approximation, like a harmonic oscillator, (jlT 



[a,a + ] = -j » 1. 

If we rewrite the Hamiltonian in terms of these operators, we find 

H = tta + a + V (a,a + ) 



(20) 



(21) 



where 



n = H + Ki + -K 2 , 



[22) 
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and the potential V is in this case equal to 

V (a + ,a) = ^K 2 (aa + a + a + ) . (23) 

In general, the spin part has the form 

H s = VLa + a + V (a + ,a) . (24) 

The calculation presented below can be tailored to deal with any polynomial V which will 
give rise to some type of Feynman rules. |l7h 

Before we end this section, we would like to point out the analogy between the spin 
problem and the harmonic oscillator problem when they are coupled to a bosonic bath. For 
a quantum oscillator, the Hamiltonian is 

n ° = m +v ® = ii + 1 2 M ^ < 25 > 

while for the spin Hamiltonian, we have 

-Hs = pl + py (26) 
S h + (27) 



2M 2 
which means 

A = Mul B = ±. (28) 
The commutation relation for the harmonic oscillator 

[q,p\=i, (29) 

and for a spin, we similarly have 



q r, 



iS z . (30) 



If 5*2 is a constant of motion, we can replace the operator S z by its average value and 
normalize the remaining components by it. Both systems, the harmonic oscillator and 
the spin, are equivalent if they are coupled the same way to the bath. Most works on 
the harmonic oscillator case involved linear coupling to the bath. Reference (|2"of ) gives an 
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exhaustive treatment of this problem. If for each bk,b\, we define two new real operators 
Xk and pk, such that 

h 



and 



Then the coupling is given by 



1 / pk 

—7= \/u k X k + I—— 

If r- . Pk 
V2 \ 



N 
k=l 



(31) 
(32) 

(33) 



The equations of motion are given by the Heisenberg equation. For the harmonic oscillator, 
we have 



iq 



q,H 



ip 



p,H 



(34) 



with similar equations for the magnetization with q — > S x and p — > S y . For the linear 
coupling, Eq.(jnHl), the correlation function of the position, 



is equal to ()2 



<?« (*) = T7 



where 



C qq (t) = -(q(t)q + qq (t)}, 



duo 



toot (—iuj) 



-2 coth — - cos [out) 



1 N 2 
v ; M ^ 9,,,, 



M ^ 2bJ k L0 2 k + z 2 ' 



k=l 



(35) 



(36) 



(37) 



and (5 is inverse temperature 1/kT with k the Boltzmann constant. For Ohmic dissipation, 
the damping kernel is without memory and hence a (z) is a constant. This is achieved 
by having a bath with spectral density linear in frequency. In this case the magnetization 
correlation function for the x— component in the high temperature limit is 



C xx (u) = 2akT — ^ 2 ; 3 • 



{ujI - uj 2 ) + (aou)' 



(3* 
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From Eq. (JSSJ), the (FMR) frequency lu is, as expected, equal to V AB. This result may 
not however be observed for a spin variable, since the interaction term has the unusual form 
of having only one component coupled to the bath, 

N 

«" = -E2^=(fc+tf)3- < 39 > 

Here it is the x— component of the magnetization which is coupled to the bath. Hence the 
coupling in Eq. seems more reasonable for a spin variable, in general. The question 
of counter-terms will not be treated here and any shift in the frequency will be absorbed. 
Reference (0) gives a detailed treatment of these terms within the models treated here. 



3 The Coherent State Representation : Equilibrium 
and Non-Equilibrium Dynamics 

Here we first review the Bosonic coherent state representation and then show how it can 
be used within the path-integral formulation of quantum mechanics. Some typical formulas 
are presented here for transition rates between two states. All the results are based on the 
standard Gaussian formula for path integrals. 

3.1 Coherent States 

Coherent states are the natural representation for semi-classical calculations. A Gaus- 
sian wave-packet for a harmonic oscillator with minimum uncertainty is a coherent state. 
They are formally defined as eigenstates of the annihilation operator (|9|) 

a | a) = a \a) (40) 

where a is a complex number. Usually the ground state is defined as the state with zero 
quanta, 

o|0) = 0. (41) 

For spin coherent states, the ground state is taken to be the state with the highest weight 
J, or zero deviation from highest weight 

\0) s =\J,J), S 2 = J(J+1). (42) 

For Bosons, the Hilbert space ( or space of all possible states) is a linear combination of all 
state vectors, {|n)}, such that 

|n> = -L=«f|0). (43) 
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These states form an orthonormal basis. Coherent states may then be written in terms of 
these states. We have 

| a) = e aa+ |0) . (44) 
Spin coherent states are defined in a similar way. We have 

1 

l + \a\ 

where the extra factor in front is due to the constraint that S 2 is constant. An important 
operator relation for a path-integral representation is the decomposition of the unit operator 
in terms of coherent projection operators 

' ,n -^ Q \a)(a\ = 1, (46) 



l«> = ^^ e l°>s> ( 45 ) 



2m 



which is used in the discretization of the path integral. (J 171 ) The coherent states form an 
over-complete basis. 

In all the calculations carried out below, we keep repeatedly using the fundamental result 
for the Gaussian integral, 

/ Udfj, (<j>i) exp I - <f>i A ij<t>j + Yl + Y1 & ai \ 

3.2 Propagators 

A typical propagator that shows up often in the calculations of the correlation functions 
is associated with a Hamiltonian that has the following general form, 

H [J, J*} = Vla + a - J(t)a + - J* (t)a (48) 
where J(t) is a time-dependent external source. The kernel of the evolution operator from 



an initial state Zi to a final state Zf is given by (jlTj) 

Kj(z f ,t f ,z i ,t i ) = (z f ,t f | Zi,U) = e Wl. (49) 

Applying the Gaussian formula, Eq. (j4T|) . with the boundary conditions 

z(ti) = z h (50) 
z(t f ) = z f , (51) 
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gives the phase S, 

s[j,r] 



-iz f e 



+ 



dt 



^ /e -*n(*/-t) + r{t)e' in ^ u) Zi 



(52) 



+i 



dt 



dt' r (t) e-™^ j (t') e (t - t') . 



O (t) is the unit step function. In using coherent states, the variables z and z are not 
necessarily related by complex conjugation. They have to be treated as independent. 

3.3 Partition Function: Euclidean formulation 

This formulation is useful to calculate thermodynamic equilibrium properties of a system. 
In this case, the partition function is obtained through a calculation of diagonal propagators 
in imaginary time, 

T = it, 0<r<{3 (53) 
In the coherent state representation, the partition function is written 



Zfs [J, J*} = / dfx(z)e- z * z ( z 



Texpl - / drH [J, J*} 



(54) 



where T is the imaginary time ordering operator. The integrals in this partition function 
are all Gaussian and hence can be easily calculated with the aid of Eq. (|4*7|) . We get 



[J, J*} = n (u) exp {(3uo + S [J, J*}} 



where 



S [J, J*] = dr dr'J*(r)D (r, r') J(r'), 
Jo Jo 

D (r, t') is the Feynman propagator of this model 

D(r,r') = (T (a+ (r) a (/))>. 
It is not difficult to find that it is given by 

D (r, t') = n (w) e-< T ' T '- f ) [ e ^9 (r - r') + e~^9 (r' 

where 



n (u) 



(55) 
(56) 

(57) 

(58) 
(59) 



is the Bose-Einstein distribution. Its Fourier transform coincides with the well known 
Matsubara propagator 

D(u n ,cj) = , u n = 2irnT, n — 0, ±1, ±2, ... . 

iuj n + to 

This last propagator is essential for any calculations that involve calculating expectation 
values of any observable, (flil) In the appendix we show how to use this method to calculate 
correlation functions. 

3.4 Real-Time Formulation: Dynamics 

The real-time formulation deals with non-equilibrium questions. In this case we can 
derive equations of motion for any observables. (fljh This is the method we adopt in the 
calculations of the correlation functions of the magnetization. For a general operator O, its 
average value at any time t is given in terms of the density matrix p, 

(O (t)> = Tr (pO (t)) . (60) 

The operator O is in the Heisenberg picture, 

0(t) = e^Oe-™*. (61) 
Therefore the average of the observable O at time t can be written in terms of that at t — 0, 

(0(t)) = Tr (pe^Oe*™) . (62) 

This latter average can be written in terms of path integrals as in the equilibrium case. First 
we define the operators K and K. The operator K is a forward propagator and is defined 
as follows, 

K [ Ji, J*] = T exp | -i jf f (h - J*a - Jia + ) d*| (63) 

= Texp|-zy (n - F*S X - FfS^j dt\ 

T is a time ordering operator. K is a backward operator and is therefore defined in terms 
of anti-ordered time operator T, 

^[J 2 ,J*] = Texpj-ijf' (ii. - J 2 *a - J 2 a + ) dt\ (64) 
= Texpj-^y (ft - F£S X - F$S^) dtj , 
10 



where F 1 and F 2 are real external fields which are coupled to the transverse components of 
the magnetization. Similar to the equilibrium case, we define a generating functional 

Z [J, J*] = Tr {p [J 3 , J*] X [J 2 , J*] K, [Ji, JJ] } . (65) 

J is now the three- vector (J 1; J 2 , J3). The density matrix p is assumed of the form 



p[ J ^ J l] = T 7 exp 
= T 7 exp 



Jo 



n 



n 



J*a 



J 3 a I dr 



(66) 



F%S y ) dr 



where T/ is now a time-ordering operator along the imaginary time axis. Hence, all correla- 
tion functions can be obtained from the coefficients of the Taylor expansion of the functional 
Z [J, J*] around J = J* = (or Fi = F 2 = F 3 = 0). For example, the average value of the 
re— component of the magnetization at time t can be found by differentiating Z with respect 
to Ff at the same time t, 



1 5Z[J,J* 



Z 5Ff(t) 



- (S x (f)> • 



(67) 



F=0 



Next we define another functional W which at equilibrium becomes the thermodynamic 
potential of the system, 

Z[F i=li2 , 3 ]=exp{iW[F]}. (68) 

The functional W, as will be seen below, is the more appropriate functional to calculate 
and expand in powers of J and J* (or F). Therefore, we have for averages and two-point 
correlation functions, 



5W 



5F? (t) 



F=0 



and 



5 2 W 



5F?(t)5F?(t') 



-i(T(S x (t)S x (t'))). 



(69) 



(70) 



F=0 



Similar expressions hold when we differentiate W with respect to the sources J and J*. They 
are related to each other by chain rule, e.g., we have 



5 5 
+ 



5F? 5Ji 5J{ 



(71) 



Depending on how we couple the bath to the spin, we use either sources to find the corre- 
sponding correlation functions. In the normal mode coupling scheme, we assume that the 
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normal modes of the spin are coupled to the normal modes of the bath. In this case, it 
is more advantageous to express everything in terms of creation and annihilation operators 
and hence use the J sources to get the correlation functions. This is what we do in section 
5. In section 6, the spin is coupled directly to the bath oscillator. Hence in this case we 
use the F sources to get the correlation functions of the spin. In the next section, we give 
an explicit expression for the functional W in terms of coherent states. 

4 Coherent State Generating Functional 

In this section, we continue working within the real-time formulation. We give the full 
expression for the generating functional in the coherent state representation and calculate 
all the associated propagators in this case. 

The generating functional Z is defined above, Eq. (J65|) . Using coherent states, for both 
the bath and spin, this trace formula can be written in terms of path integrals over spin 
variables and bath variables, 

Z [J, J*] = J d/i(a 1 )dii(a 2 )dii(a 3 )dii((p 1 )dii(ip2)dii(ip 3 ) (a u ip x \p [ J 3 , J 3 *] | a 2 , <p 2 ) 

x (a 2 ,y?2 |/C [J 2 , J 2 *]| a 3 ,v?3> (a 3 ,<Pa \K [Ji, JJ]| a lt </?i) (72) 




The aj=i,2,3 represent states of the spin system, while the ^=1,2,3 represent the bath 
states. This integral can be formally written as a path integral along the path in FigQ 
with periodic boundary conditions similar to the equilibrium partition function calculations. 
This functional can be calculated exactly only in few cases in particular if the Hamiltonian 
is quadratic. Higher order terms can be accounted for only approximately. This is best 
done through a graphical procedure such as the Feynman diagram technique. Here we have 
a quadratic Hamiltonian and hence we can solve for Z, however we will mention briefly what 
happens in the general case. 

In our case, the bath degrees of freedom can be integrated out exactly and we can derive 
an exact effective action for the spin degrees of freedom. In the general case, the effective 
action can be derived perturbatively. From it, we will calculate the correlation functions of 
S. We find 
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Z[J,J*] = J dfi(ai) J dfj,(a2) J dn(otz) exp {— |«i| 2 — |«2| 2 — |«3| 2 } (73) 

/■«3 r<*2 pcil f ^ I 

x / dfi(zi) / dfi(z 2 ) / dfifa) exp < ^ J» [z h z h J i: J*) > T (Z, Z) 

</ai J03 i«2 ^ j = j J 

where JF (Z, Z) is the Feynman- Vernon functional for the spin-bath system. It is given by 



ln^(Z,Z) = J dt J dt ' 



(74) 



where the three-vector Z is related to the three branches of the curve C, Fig. Q 



£1 

Z3 



(75) 



The time integrations are defined based on the path C: 



<t,f <t 



/• 



t,t' e CW.CH 



t» < t, t' < u - ip, t,t'eC {0) . 

The Feynman- Vernon term is the only term which is dependent on the bath parameters. 
The functions (t,f), nine in total, are propagators associated with the bath oscillators. 
Hence they can easily be calculated using Eq. (152)) since the oscillator part of the Hamiltonian 
is quadratic and the spin can be considered as the external field. The indices i,j = 1,2,3 
relate to the branches C<+>, C<® of C. They are (Q) 



6 (t - t>) G k 21 + 6(t- t) G k 2 

n(uj k )e^ k e iUki - t+iT) 



n (u k ) e - <w *<*-*') 
0(t-tf) G k 2 + 9{t' -t)G k 21 
GUt-r) 



io;fe(t+ir) 



n (u k ) e 



-U3 k \T—T 



n[uj k )e 
G k 13 (t,r) 

1 \6 (r - t') e"^ 



+ 0(t' - t) 2 



(76) 
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These "path-coupling" functions show that if the initial time i, is taken to be in the infinite 
past, ti — > — oo, the branch decouples from the other two branches, (flih This is 
the case where transient effects have died out. In the rest of this paper, these transient 
effects will be neglected and we will concentrate only on the real-time paths and C^~\ 
We take account of the third branch through the assumption that initially the system is in 
equilibrium. For a general potential V, the generating functional can be written in terms 
of that of a free system, 7i = fla^a, interacting with the bath, 

d 8 



Z [J, J] 



exp 



-i I dtV 

c 



%SB [J, J] 



(77) 



dJi{ty dJi{t). 

7isb is therefore the generating functional of a particle interacting with the bath and there 
is no external potential. This latter formula is valid in the general case and is the start of 
any perturbative calculations. The action along the real-time trajectories is given by 



ill 



dt 



along the path and by 



i/° = -i dt 



Z\Z\ — Z\Z\ 

2i 



Z2Z2 — Z2Z2 



flZiZi 



(7S 



2i 



- Qz 2 z 2 



(79) 



Z{U) 



-pH(U) 



along the path \ Fig. [TJ At U — > —00, the system is at equilibrium, then we can assume 
that the initial density matrix is thermal, with J^{ti 

p{ti) = 

Then, we observe that 

(an oo)| a 2 ) 



= 0. Therefore we write that 

U -> -00 (80) 



dfi(z 3 )e u °^^(z 3 ,z 3 ) 



«i 



where J° has the same expression as Jj 1 but with t — > it. Hence in this case, the initial 
density matrix element is just another overall factor in the generating functional Z, 



Z [J, J ] 



d[i(a 3 ) / dfjL(ai)dfj,(a 2 ) (aii \p(— 00) | a 2 ) 



(82) 



x 



x 



03 



d/i(zi) 



«3 



dn(z 2 ) exp 











,. + dtV \dJ(ty dJ(t) 



exp [ Zl , zt] + ill [Z2,z 2 ]+i J dt (J • Z + J • Z) W (Z, Z) , 



14 



with Z = (zi, —z 2 ) and the measure is defined by 

d]I (a) = dji (a) e~' Q ' . 
Therefore we define a new generating functional Z 



(83) 



Z [3,3] 



dfi(as) / d~jl(oii 



dn(a 2 ) (cti |p(-oo)| a 2 ) x Z [J, J] 



We can now adopt a different notation that takes into account the path C implicitly by 
defining a scalar product and combine the different components into a single vector. The 
generating function becomes 



Z [J,J] 



dfi(zi) / dfi(z 2 ) exp 



dtV 



d 



d 



d3(tYd3(t) 

xexp (Z,Z)+i y rft (Z- J + J-Z)|j !r (Z,Z) , 

where now the vector Z is defined 



^2 



and the free action is 



with 



/ = E-^ 



1=1,3 



1 
-1 



The complex scalar product is now defined by 



Z J = a lJ z]Jj. 



(84) 

(85) 
(86) 

(87) 
(88) 



This notation makes it possible to take into account the closedness of the real-time path 
by just taking one branch of the curve C and doubling the components of the dynamical 
variables. The matrix plays the role of a metric. 
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Now we turn to some properties satisfied by the functions G k j. These properties are better 
displayed in Fourier space. The Fourier space representation of the Feynman propagator is 
given by 

G k u (uj) = J dt e iujt (l+n (w fc )) 9 (t) e~ iuJkt + J dt e tuJt n (u k ) 9 (-f) e - *"** (89) 
= (l + n(u k )) l — -n(u k ) -, (90) 

UJ — U3 k + 26 UJ — U) k — l€ 

where e — > + . V stands for the principal part of the integral. For the anti-time ordered 
propagator, we have 

G k 22 {uo) = ! dte iut n{uo k )Q{t)e- lWkt + [ dte iujt (1 + n{u) k ))Q {-t) e~ iulkt (91) 



= n{u k ) —-(l + n(u k )) -. (92) 

UJ — U) k + It UJ — UJ k — l€ 

For the other remaining Green functions, we have for positive uj 

G k 2 (uj) = 2-nn (uj k ) 5 (uj - uj k ) , (93) 

and 

G k 21 (uj) =2n(l + n(uj k ))6(uj-uj k ). (94) 
These Green functions are not all independent. We first observe that 

G k u (uj) + G k 22 (uj) = G k 2 (uj) + G 21 (uj) . (95) 

This is an immediate result of their definition. Moreover, the term on the l.h.s. is easily 
seen to be a symmetric sum of the product of two operators. Now it is not difficult to show 
from the above expressions of the Green functions that we have 

G k u H + G k 22 (uj) = (1 + 2n(u;)) [G 21 (uj) - G k 2 (uj)] . (96) 

The last factor on the r.h.s. is an anti-symmetric sum of two operators. Equation (J96|) is 



a statement of some form of the fluctuation dissipation theorem. (|2l[). These relations will 
be used in subsequent sections to calculate the correlation functions. In equilibrium, when 
the distribution functions are the Bose-Enstein functions 

1 + 2n(uj) = coth — (97) 
2 
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and equation (pj)) gives the usual from of the fluctuation-dissipation theorem. 



5 Normal Mode Analysis 



In this section we follow Lyberatos, Berkov and Chantrell (|g) and couple the normal 
modes of the spin system (also called collective field) to the harmonic oscillators of the 
bath. This method has also been recently used by Safonov and Bertram to calculate 
correlation functions of the magnetization in thin films. (0) In this section, we show how 
their correlation functions for the collective field follow from the microscopic model treated 
here. The results in this section will be used in the next section to find the correlation 
functions of the magnetization in the general case. 

First we find the collective degrees of freedom c and c* from the magnetization: 



a = uc + vc* (98) 
at 

where u and v are real. We require that 



a) = vc + uc* 



This implies that 

Therefore, we can write for some 9, 



If we set, 



4 

we find that the coefficients of the transformation Eq. (l98|) 



[c,ct]=l. (99) 
u 2 -v 2 = 1. (100) 



u = cosh# (101) 
v = sinh 9. 



ul = n 2 -^, (102) 



u = ^~2^T' (103) 



Q — uoq 
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In this collective coordinates, the spin Hamiltonian becomes diagonal, 

Now the interaction term between the bath and the spin is taken of the form 

Hsii = -5](7A + #). 



(104) 



(105) 



The generating functional for this system is then calculated in Fourier space with the help 
of the Green functions stated in the last section. We find that 



Z 



[J, J] = j dfi (z u z ± ) J dfi (z 2 , z 2 ) exp j- J^- [Z-^.Z + J-Z + J.Z]| 



det.4 



exp 



2n 



where Z = (z\, z 2 ) and J = (Ji, J 2 ) are two-component vectors. The matrix A is 



A 



13 



i (uj -u) + rin M n i2 (u) 

n 2 i (u) -i (uj — u) + n 22 (00) 



(106) 



(107) 



The Uij terms are due to the interaction of the system with the bath. They depend on the 
density of states of the bath A(o>) and the coupling constants. For a general bath, the Pin 
element is given 



IlnH = X>*| 2 GjiM 



J ~ 



VTA (U k ) |7 M 



1+n (u k ) 



n(u k ) 



(108) 



(109) 



uj — uj k + xe u — uj k — xe 

The remaining matrix elements are calculated similarly. If now, we assume that the bath 
parameters satisfy the condition 

7rAK)| 7 K)| 2 = a (w fc >0), (110) 

where a is a constant, we find that the interaction with the bath induces the following 
coefficients, 

U n (co) = a(l + 2n(w)) (111) 

U 22 (cu) = n n (w) (112) 

II12H = 2cm (w) (113) 

n 2 i(w) = 2a(l + n(w)). (114) 
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The correlation functions of the magnetization are related to the inverse elements of the 
matrix A. A calculation of the inverse matrix gives 



1/1 V 



i (uj — a;) + a (1 + 2n(u>)) —2an (uj) 

-2a(l + n(oj)) -i (oj - u) + a (1 + 2n(uj)) 



where T> is the determinant 



D = dct^l= (uoq - uj) 2 + a 2 



(115) 



(116) 



After solving for the generating functional in terms of the external sources, we can expand 
it around the point J = J* = to get all the required correlation functions. In this case, 
we can solve for Z [J, J] exactly since the full Hamiltonian is quadratic. The correlation 
functions are now found by differentiations with respect to J and/or J. For t > t', we have 
for the correlation function of the collective operator c 



J=J=0 



where 



W 



[ J, J] —i J 



duj 



117) 



(118) 



All correlation functions of three operators or more are zero since the Hamiltonian is quadratic. 
The above correlation function is therefore related to the matrix element, A^ ■ At high 
temperature, i.e., (3 — > (uj « kT) 



A^(uj) = 



2an (u>) 



(uj - uj) + a 2 
-2akT 1 



UJ 



(ljq - uj) + a 2 



In real time, we have 



W 



[J,j]=ijdtj dt' Ji (t) A' 1 (t - t') Jj (t') , 



(119) 
(120) 

(121) 



hence, the two-point correlation function of the field c associated with the 1 — 2 component 
of the path C is 

27r (^ — uj) 2 + a 2 



(J(t)c(t')) = 2a 



(122) 
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Similarly, we can find the corresponding 2 — 1 component of the correlation function, 

From these last two correlation function, we get the classical correlation function of c, 

\ (cHtW) + c( t ')ct(t)> = a / 1 + 2 " ( "' - . (124) 

2 y 27T (a; - a;) + a 2 

Now we observe that if we take the limits a — > and t' — > t, we recover the expectation 
value of the occupation number n = cc^ , 

(n{u)) = l + n{u). (125) 

To get this limit, we have used the fact that 

-^-^ — - — > 5 (u) as a — > 0. 

it uj z + a z 

We also note that in this limit, we have 

(c f c + cc f ) = 1 + 2n (to) , (126) 

and 

<[c(t),ct(t)]) = l. (127) 

This shows that the commutations relations and the FDT are satisfied at all times t. 

Since we are close to equilibrium and a is small, the power spectrum will be picked near 
uj = ujq. Therefore, this is also equivalent to having a Langevin equation with random forces 
F such that (19) 

c = -(iu + a)c + F(t) (128) 

and 

({Ft (*')}) = 2akT5(t-f). (129) 

Within this approximation, we recover the correlation functions of Safonov-Bertram, Eq.(2.14) 
in Ref. (0). 

{c\t)c 0) = / -TT-—. rj— -• (130) 

' J Ztt uj (cuq - uj) + a 2 

To get the correlation functions of the magnetization S, we first use the linear transfor- 
mation, Eq. (J98j) . to write S x and S y in terms of the collective operators c and c'. This 
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result will be deduced in the next section from the exact treatment of the general asymmetric 
case and without recourse to the rotating wave approximation as was done in Ref. (0). 

6 White and Colored Noise: LLG and Other Solutions 



One could compute the correlation functions for the original variables from the cc^ correlation 
functions, by inverting the Bogoliubov transformation ([98)1 . However, in this section, we will 
repeat the computation directly in terms of original magnetization component variables, 

S x = i(a + at), (131) 

S v = ^(a-O • (132) 

Writing the generating functional Z in terms of them is trivial, but performing the Gaussian 
integration is more complicated, since we will have to invert a 4 x 4 matrix. We couple the 
magnetization to an external time dependent magnetic field F 

H e — — F ■ S. (133) 

As we have seen in the previous section, in the generating functional approach, we double the 
components of S and that implies doubling of the external field F. Therefore the interaction 
term becomes 

# e = -F 1 -S 1 + F 2 -S 2 . (134) 

The results for this Hamiltonian can be derived from those already found in the previous 
section. The bath contributes a term of the following form to the effective action 

S eff (S 1; S 2 ) = J ^ | 7fc | 2 Gu M K, + Sl,y + iSi, x Si, y - iSi, y Si, x ] (135) 

- j — \lk\ 2 G\ 2 (u) [S% iX + S^y + iS 2iX S 2 ,y - iS2,yS 2 ,x] 

/doj 
2~~ 1 7*! I G 21 (id) [S2,x.Sl,x + S2,ySl : y + iS 2 ,xSl,y ~ iS 2 ,ySl >x \ 

+ j ^ ^fc^ ( U ) [Sl,xS 2 ,X + Sl t yS 2> y + iSl >x S 2j y ~ iSl^yS^^ ^ 
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where the bar denotes the complex conjugate integration variables and Si (S 2 ) is the com- 
ponent along the path (C- - - 1 ), Fig. ^ Next, we define two new vectors S and D, 



Similarly, we define 



S 
D 

F d 
F„ 



i(s 1 + J 

Si — S2. 



F,-] 

i(F, 
2 v 



Finally, we make another definition. We define four-vectors U and F 

U = [S X , Sy, D x , Dy) , 

F = (F d ,F a ), 



(136) 
(137) 

(138) 
(139) 



(140) 
(141) 



and write the generating functional in terms of these four-vectors along the path C^' only. 
Since U (t) is real, then we have 

U(uj) = U(-uj), 

and hence we should constrain the fourier integration to positive frequencies only. The 
bath-independent part of the Hamiltonian then gives the following contribution to the phase 
of Z, 



Hi - ih 



duj 

7T 



u) Uj {to) , 



(142) 



where the matrix is, in Fourier space, 



(0) 



iA -uj 

u iB 

iA -co 

uj iB 



(143) 



Again, it is the inverse of the full matrix A = A^ ' + A mt , that is needed to determine 
the correlation functions of the magnetization where A mt is the part that is due to the 
interaction with the bath. The determinant of A determines the natural frequency of the 
system and any broadening due to interactions. The determinant of the free part is 



2\ 2 



V = (uj 2 - uj 2 ) 



(144) 
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where 

ul = AB, (145) 

is the ferromagnetic resonance (FMR) frequency of the system. The calculation of the 
matrix A is done along the same lines as in the normal mode solution. 

To recover dissipative behavior in the spin sub-system, we take the continuum limit in 
the number of oscillator modes. This limit guarantees that the probability of acquiring back 
any energy lost to the bath is zero. Because of the interaction with the bath, we expect 
that there will be a shift in the energy of the spin system accompanied by dissipation. 

An explicit computation shows that in the continuum limit, i.e. converting the sum over 
k in an integral over the frequencies involving the density of states X(uo k ) 

X(u k ) = ^- , (146) 

the correlation functions can be expressed in terms of the functions L r and Li : 

L r (w) = -i r ^ ttA ( Wjfc ) | 7 M\ 2 (G k u - G k 22 ) (u k ) , (147) 

Jo TT 

POO 1 

U{u) = 2/ — ttA (uj k ) |7 (uj k )\ 2 (G k 12 - G k 2l ) (w fc ) . (148) 
Using the definitions of the Green functions, we find 

roc j i 

L r (oo) = 2V — vrAK)| 7 (^)| 2 , (149) 

Jo 71 uj -u k 

Li(u) = -2n\(cu) \-i (uj)\ 2 6(uj) . (150) 

Counter-terms are needed to cancel ultraviolet divergences in L r . For simplicity, we will 
assume that this is done via suitable subtractions. The effect of L r is a redefinition of the 
given coefficients A and B. This redefinition in principle changes the oscillation frequency. 
However, for a passive path, one neglects L r (u) and thus the frequency shift. In this approx- 
imation the coefficients A and B are kept unnormalized and all the physics is contained in 
Li(u). There is a subtlety here, since the expression (|15()jl is not antisymmetric, whereas it 
has to be antisymmetric due to general properties of correlation functions (see appendix for 
a discussion). Therefore Li(u) has to be antisymetrized. By noticing that ^(w)] 2 is even in 
uj and extending X(u) to negative u < with a negative sign, the final result can be written 
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in the form 

iA-A{u) -uo 

uo iB-A{u) 

iA + A(u) -co n\(uo)\-f(uo)\ 2 (l + 2n(uo)) 

u iB + A{uo) it\{uo)\~i(uo)\ 2 (l + 2n(u)) 

(151) 

where A(u>) is the odd function 

A M = =hM+hW = * X{uMu) r. (152) 

The determinant of this matrix is given by 

detA = V(u) = [uo 2 -uo 2 - A(wf] 2 + [A(uo)(A + B)] 2 , (153) 

We also observe that for the functional integral to converge, we must have 

n\(uo) |7H| 2 (1 +2n(u)) > 0. 

This requires that the function Li (u) when extended to negative frequencies be an odd 
function which is consistent with the statement before Eq. (jl51|) . 

To calculate the correlation functions, we need first to calculate the inverse matrix of 
A. The cofactors needed for the correlation functions of the different components of the 



magnetization are for small couplings to the bath 

Cl i = -ttA (u) |7 Ml 2 (1 + 2n(u)) [B 2 + lu 2 + A (lu) 2 ] , (154) 

C12 = | 7 H| 2 (l + 2nH) [(A + B)u], (155) 

c 22 = -vrA(cu) | 7 (^)| 2 (l + 2n(w)) [A 2 + u 2 + A(u) 2 ] , (156) 

C12 = -c 21 = i7r\{uj)\ 1 {uj)\ 2 {l + 2n{uj)){{A + B)uj]. (157) 



As will be seen below, the en (C22) co-factor of the matrix A is associated with the 
xx (yy) —component of the magnetization while ci 2 is related to the xy— component. 

Correlation Functions 

Next, we use these cofactors to calculate the correlation functions of the magnetization. 
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For a general operator O, the average of the anti- commutator {0(t), 0(t')} is found by 
differentiation of W [F a , F d ] with respect to F d , 

\ { 0(t)0 V ) + 0(f)0(t)) = -i ™^^ (158) 

Applying this procedure to the components of the magnetization, we find that for the 
x— component 

\{S x {t)SM + SMSS))= I ^cos^^M. (159) 

l J ATT V [00) 

From Eq. ()154|) . we then obtain 

C xx (t) = I ^cosM)[l + 2nM]7rAM| 7 M| 2 ^! + ^l±^! , (160 ) 

Now, we show how for different choices of the function Lj (to), we can recover the LLG result 
and other oscillator-type correlation functions. 
6.0.1 Case 1 : LLG 

If we assume that the bath is defined such that 

7rA (oo) |7 (uj)\ 2 = au, (161) 

i.e., Li (00) is odd. Then in the limit of high temperature, f3 — > 0, the correlation function 
for the xx-component takes the simple form 



^ cos(^t) 



(162) 

\(l + a 2 )u*-uff + [atv (A + B)] 

This is the result that coincides with that derived from LLG. (0; 0) This case also corre- 
sponds to a white noise solution, (flih Moreover, we observe that the condition on the bath 
that gives LLG is similar to the one that gave the harmonic oscillator solution, Eq. (j38|) . In 
both cases the spectral density is linear in frequency. (20) 
6.0.2 Case 2 : Coherent Oscillator 

This case is similar to the normal mode result. We call it coherent oscillator because 
this case gives correlation functions similar to those of the collective operator c in the normal 
mode analysis. Here we choose Li(uj) such that 

n\(uj) |7(cj)| 2 = a, to > 0. (163) 
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This is the same choice as in the previous section. For (3 
for the ^-component of the correlation functions 



0, we get the following expression 



C xx {t) = 2akT 



du ( 1 

— coslcut) — 



B 2 + u 2 + a 2 



Co> z — LU, 



If 



+ 2a 



2 (J 1 + w, 



+ a 4 



(164) 



The normal mode result is easily seen to follow by setting A = B in the correlation functions 
and (1/uj) by (1/uj ) . This result has been obtained before by Safonov and Bertram. (p) 
However without this latter approximation, this model corresponds to a case of colored 
noise. (Jilt ) As ui — > 0, the integral diverges. Therefore at small u, the approximation in 
Eq. (|163|) is not applicable: Li(u) — Lj(— cu) cannot be a constant but, for consistency with 
antisymmetry and analyticity, must vanish with u at u — > 0. 

7 Conclusion 

Starting from simple quantum models which only differ in how the spin couples to the 
bath, we have been able to derive variant correlation functions for the magnetization close 
to equilibrium. We have limited ourselves only to linear-type couplings. Depending on the 
coupling and the density of states of the bath, we showed how to obtain different types of 
correlation functions including the classical LLG result. First, we showed that the typical 
harmonic oscillator correlation functions are recovered only if the S x component of the mag- 
netization is coupled to the bath oscillators. Next, we coupled the normal modes of the 
spin to those of the bath and this allowed us to get the correlation functions of the collective 
field which is the starting point of the work of Safonov and Bertram. We were also able 
to use this special coupling to derive a more general type of correlation functions without 
recourse to any approximations such as the rotating wave approximation. These correlation 
functions are for a general linear coupling between the bath and the magnetization depend 
on the coupling constant 7 and the density of states of the bath system, A (a;). For the 
SVSV-correlation function we find, 



C xx {t) = I [1 + 2n(u)] A (u) cos(ut) 

Z7T 



B 2 + A (cj) 



(u 2 - uj 2 - A {uo ff + [A (w) (A + B)\ 



Similarly for the S'^S'y-correlation function, we have 

/d 
— [1 + 2n(u)] A (u) cos(cut) 
2% 



uj 2 + A 2 + A M 



(u 2 - uj 2 - A {uoff + [A (u) (A + B)f 



(165) 



(166) 
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and finally for the S x S y -coTTe\aXion function, the correlation function is 



Cxy(t) = / — [1 + 2n(w)] A (a;) sin(wt) 




a; (A + 5) 



(w 2 - a; 2 - A (^) 2 ) 2 + [A (w) (A + B)] 2 



where 



A (a;) = 7rA (u) (7 (co>) | 2 . 



(168) 



The LLG solution was obtained for a special type of density of states and coupling to the 
bath. The same condition was also obtained in Ref. (flih where in addition we were able 
to show that this choice gives the white noise character in the stochastic formulation . The 
normal mode solutions are however in general with memory. The assumption that damping 
is constant close to the FMR frequency makes the equations of motion Markovian. (|22T) The 
damping in the cases treated here is independent of the symmetries of the Hamiltonian spin 
system, the reason being that the dissipation kernel only depends on the coupling with the 
bath and the bath properties, but not on the spin Hamiltonian. For couplings other than 
linear, the damping is expected to depend on the symmetries of the full Hamiltonian, but, 
again, not at the leading order in perturbation theory. The point is that for non-linear cou- 
pling, the effective Hamiltonian and therefore the correlation functions have to be computed 
perturbatively in terms of Feynman diagrams; in particular the spin propagator will enter in 
higher order computations and since the spin propagator depends on the symmetry of the 
spin Hamiltonian, which could be isotropic (A = B) or not (A 7^ B), we will have different 
results for A = B and A ^ B. This does not happen at leading order in the non-linear case, 
and does not happen at all orders in the linear case, where the exact result is available. 

One last word about the non-equilibrium machinery used here to derive the above corre- 
lation functions: this choice of method allows us to go beyond the equilibrium formulation, 
and in particular to show that a generalized fluctuation-dissipation theorem (JSfij) holds true 
even if the distribution functions are not exactly the Bose-Einstein ones. In principle, an 
analysis of this system when the distribution functions present strong differences from the 
thermal one is also possible in the general formalism we discussed here. However, such a 
strongly out of equilibrium analysis would require further study and it is outside the scope 
of the present paper. 
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Appendix 



Here we briefly show how the correlation functions derived in the main text can be derived 
using the equilibrium imaginary time formalism |l7h and without going to the coherent state 
representation. This is an equilibrium consistency check for our non-equilibrium computa- 
tion. 

The basic idea in the equilibrium computation is to invoke the fluctuation-dissipation 
theorem (which is not assumed in the non-equilibrium computation) in order to derive the 
fluctuations from the dissipation, i.e. from the spectral density. The fluctuation-dissipation 
theorem says that in equilibrium the symmetric correlation functions can be written in terms 
of the spectral densities p l ^{oj) as 

< {S\t), S j (0)} >= J duj e iuJt coth ^ p ij (oj) , (169) 

where i and j denote the indices x and y respectively. From this definition it is immediate 
to see that the spectral densities must satisfy the relationships 

p%0* = p*(w), p^(uj) = -^(-uj). (170) 

In particular, p xx and p yy are real and antisymmetric: 

= p» = -p*(-aO- (171) 

Thus, one can extract the spectral densities p ll (uj) from the spectral representation of the 
retarded self-energy, 

D%{uo) = f duo' ptJ(u/) (172) 



by taking the imaginary part: 



p» = -- Im D%(u;) . (173) 

7T 



Moreover, due to the behavior of the theory under time reflections t — > — t, 

< {S x (t) , 5^(0)} >= - < {S x (-t) , S x (0)} > (174) 

we have that p xy (uj) and p yx (u) are imaginary and symmetric: 

p^)* = -p^H, p i i(u)=p i i(-u), i^j (175) 
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As a consequence, the p lJ (u;), i 7^ j spectral densities can be extracted from the real part of 
the retarded self-energy: 



7T 



(176) 



In order to compute the retarded propagators, one has to compute the Euclidean effective 
action obtained by integrating out the bath degrees of freedom in the Euclidean functional 
integral 



where 
with 



[db* k dbk] exp 



/ drL E (S\b k ,bl 
Jo 



L*(S\ b k , b* k ) = L% (S*) + L%(b k , b* k ) + LUS\ h, b* k ), 



Lf(^) = S x id T S v + ~A(S x ) 2 + ^B(S v ) 2 , 



L E R {h,bl) 



J2b* k (d T -u k )b kl 



L E SR (S\b k ,bl) = J2 b l^ S - + S ^- 



(177) 
(178) 

(179) 
(180) 

(181) 



Since the integration on b k and b* k is Gaussian, S e ff can be computed exactly and is quadratic 
in the spin fields: 



dr I dr' - S\t) D 



T 



t') S j (r'). 



(182) 



The Euclidean propagator is easily obtained in the Matsubara formulation by inverting the 
2x2 matrix 

n -i/ \ ( A u n \ fU E (u n ) \ , . 

where the first matrix is the inverse free propagator and the second matrix is the self-energy 
matrix (to be compared with the real time result Eq. ([108)) ) 



U E (u n ) = 2j2hk\ 2 G h E M 



;i84) 



and where u n = 2-KnT, n = 0, ±1, ±2, . . . are the Matsubara frequencies. The inversion is 
trivial. The retarted propagator can be obtained with an analytic continuation u n — > iuo: 



1 (B + U R {u) -zu 
DM V «*> A + U r (uj] 



185) 
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where T>(u) is the determinant 

V(u) = oo 2 - u 2 + (A + B)U R (oo) + U 2 r (uj) 
In order to compute real and imaginary parts, it is convenient to split 

U R (u) = Re IL(u) + i Im 1%) 
and to introduce the real quantities 

A(u) = A + Re II(w), = B + Re II(w) 

Then the inverse determinant reads 

1 



(186) 
(187) 



Loi-u 2 -\Il R \ 2 + (A + B)U 



R 



UJ 2 — u 2 — 




2 + (A + 5)( Re n-i Im IT) 







with 



|P(a;)| 2 = [uol -co 2 - \U R (u)\ 2 + Re U(u)(A + B)} 2 + {A + B) 2 (lm U(u)f 



(189) 



(190) 



The function Re II (a>) and Im IT (a;) are related to the previously defined functions L r (ui), 
Li(uj) and A (a;). In particular 

Im n(a>) = — 27r |7 fc | 2 5(c<j — uok) = —tiX(lj)\j(lj)\ 2 



-A(uj) . 



(191) 



Notice that the continuum limit has been taken by ensuring the antisymmetry of Im U(uj). 
The p n (u>) spectral densities are obtained by taking the imaginary part of the full retarded 
propagators D R (cu): 



1 



7r|X>| 
1 



7T 



LUn — LU 



LUn — LU 



\U R \ 2 - B(A + B) 



\U R \ 2 -A(A + B) 



A 



A 



(192) 



(193) 



whereas the spectral densities, p l ^{oo) (i ^ j), are obtained by taking the real part of D R (u): 



u(A + B A) 



(194) 



These results coincide with equations ()165l 11661 I167|) . Therefore there is full consistency 
between the real time and the imaginary time formalism. 

The LLG limit for small damping is recovered when Re II — ► const, Im II — > a u, and 
the coherent oscillator is recovered in the region \u\ ~ u)q when Re II — > const, Im II — > 
a sgncu. The term Re IT is usually set to zero after being absorbed in the definition of the 
FMR frequency. 
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Figure 1: Complex time-path for the generating functional 
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